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We propose a general setting for a universal representation of the quantum structure on which 
quantum information stands, whose dynamical evolution (information manipulation) is based on 
angular momentum recoupling theory. Such scheme complies with the notion of 'quantum simulator' 
in the sense of Feynman, and is shown to be related with the topological quantum field theoretical 
approach to quantum computation. 
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INTRODUCTION 



In the past few years there has been a tumultuous activity aimed to proposing novel conceptual schemes of in- 
^ , terpretation of quantum computation. Curiously enough, most of them are based on topological notions. Among 
these, anyonic quantum computation [Q, fermionic quantum computation j^, localized modular functor quantum 
' field computation holonomic quantum computation |^ have mostly attracted attention. Such models appear to 
■ be simply different realizations of a unique conceptual scheme that incorporates all of them as particular instances. 
In this note we aim to claiming that such schemes may all be identified with a model of quantum simulator (in the 
sense of Feynman Q) based on (re)coupling theory of SU{2) angular momenta (see, e.g. 0| and references therein). 
] The paper raises a number of issues without entering in too many technical details but rather trying to establish the 
(~| • guiding philosophy, and is therefore foundational. 

Oh! The scheme automatically incorporates all the essential features that make quantum information encoding so much 
■^Ij • more efficient than classical: it is fully discrete (including its time-like variable); it deals with inherently entangled 
^ I states, and thus incorporates all achievable complexity in its set-up, which is naturally endowed with a tensor product 
^ ' structure; it allows for generic encoding patterns. 

In ref.|^ Feynman lists a minimal set of requirements as essential for the proper characterization of an efficient 
J" I quantum simulator: i) locality of interactions; ii) number of 'computer' elements proportional to (a function at most 
polynomial of) the space-time volume of the physical system; iii) time discreteness (time is 'simulated' in the computer 
by computational steps). Our argument is based on the fact that all such basic features are typical of spin networks. 



^ [ By 'spin networks' we mean here - contrary to what happens in solid state physics, but somewhat in the spirit 
- - ■ of combinatorial approach to space-time representation ^ - graphs the node and edge sets of which can be labelled 
by quantum numbers associated with SU (2) irreducible representations and by SU (2) recoupling coefficients, respec- 
tively. Spin networks can thus be thought of as an ideal candidate conceptual framework for dealing with tensorial 
transformations and topological effects in groups of observables. The idea is to exploit to their full extent the discrete- 
ness hypotheses ii) and iii), by modelling the computational space in terms of a set of combinatorial and topological 
rules that mimic space-time features in a way that automatically includes quantum mechanics. 



THE SIMULATOR COMPUTATIONAL SPACE 



We begin by defining the structural setting of a universal quantum simulator 9Jl which satisfies all axioms proper 
to the quantum Turing Machine 9Jt, whose computational space is identified with a spin network, can encode 
information, undergo unitary transformations, and simulate any finite quantum system completely described by 
eigenstates of SU{2) angular momentum operators. In the sequel we shall indicate how such a scheme can be 
extended to mixed states. 
Coding Information 
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The machine 9Jl building blocks are an ordered collection of n+1 mutually commuting angular momentum operators 
{Je \ £ = (for example associated with a set oin + 1 kinematically independent particles), with eigenvalues 

parametrized by ji, . . . , with j( ~ 0, i, 1, |, . . . . Such operators are assumed to sum to a sharp total angular 

momentum J with projection J^, whose quantum numbers are, respectively, j and m, ~j < m < j in integer steps. 

For any given pair {n,j), all possible binary coupling schemes of the n + 1 angular momenta ji together with the 
the quantum numbers /ci, . . . , fc„_i corresponding to the n — 1 intermediate angular momenta, provide the 'alphabet' 
in which quantum information is encoded. The resulting Hilbert spaces 7i;^(fci, . . . , kn-i), each (2j + l)-dimensional, 
are spanned by complete orthonormal sets of the form 

■ ■ • , jn+i ; fci,...,fc„_i; j,m) EE |b}} . (1) 

Such states can be pictorially represented by rooted labelled binary trees in which each node corresponds to an 
angular momentum quantum number: the root of the tree to j, the internal nodes to the intermediate fci, . . . , 
and the terminal nodes to ji, . . . ,jn+i- An equivalent representation is the binary bracketings notation proposed in 
0. Fig.0 shows an example of these kinds of representation where, fixed an ordering ji, j2, • ■ • ,jn+i and given a 
common j, there exists a correspondence between states given in (j^) and the (equivalent) combinatorial structures 
represented by binary bracketings and labelled binary trees (cfr. e.g. [^). 

Fig. 1 

Notice moreover that to each of these assignments (states (|l|)), there corresponds a unique non- associative structure 
over the tensor product W^® - ■ = span{|jimi)®- • - (g) |j„+im„+i)}, which is manifestly intrinsically entangled. 

A code, in this picture, is a sort of generalized quantum Godel number in the base associated with a field made of 
the ordered labels of the intermediate angular momenta, which range over a finite domain {e.g. for Ji + J2 = J12, 
lii — J2I < ji2 < ji + J2), as weU as the coupling brackets. Such coding spans (and defines) the space of all possible 
computational states. 

Operations as Unitary Transformations 

In the structure described above, any quantum operation is implemented by some transformation connecting pairs 
of binary coupled states, namely by the so called 'recoupling coefficients', or inj symbols 0. 
Indeed the 3nj symbols are unitary probability amplitudes 

= O'l, • ■ ■ ijn+i ; k[,.. . ,fc^_i ; j,m 1 ji, . . . , j„+i ; fci, . . . , A:„_i ; j,m) , 

\^3nj\^ representing the probability that the system, once prepared in state . . . , j„+i; fci, . . . , fc„_i; j, m), is mea- 
sured in state . . . , j„+i; k[,.. . ,k'^_^;j,m). 

Notice that the recoupling coefficients can be interpreted as reduced matrix elements since the total magnetic 
quantum number m can be neglected in view of the Wigner-Eckart theorem. Since they give the elements of the 
transfer matrices connecting any pair of states, the symbols actually provide the (matrix) analog of the transition 
function of the quantum Turing Machine [||. 

Moreover as any transfer from a state to another - states being in one-to-one correspondence with the vertices of 
a suitable graph (S„ as we shall see below - can be thought of as a (discrete) path integral, 3nj symbols implicitly 
define an 'action' (and hence an associated hamiltonian operator). 

We shall show in the next section that programming 971 consists just in selecting which transformations do perform 
the desired computation. 
Computational Space of the Machine 

The computational space associated with 071 is a graph C5„ whose vertices are identified (i.e. are in one-to-one 
correspondence) with the system pure angular momentum eigenstates defined in (^. 

The Racah transform TZ, together with the phase transform 

^ : |...((a6),c)^...)^|...(a(6c)J^...), (2) 
$ : \...{ab)...)^\...{ba)...) , (3) 

exhaust all possible transformations between pairs of binary coupling schemes for any n. We shall refer to this 
statement as Biedenharn-Louck theorem (topic 12 in 0). Interpreted as transformations on binary trees, TZ and 
represented pictorially in Fig.^ and Fig.|[ are referred to as rotations and twists, respectively. 
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Fig. 2 
Fig. 3 

The coding proposed above requires both types of operations, (^) and (||), and the corresponding graph is the full 
twist-rotation graph ©„, the vertices of which are to be associated with the computational states of the Machine and 
the bonds with either Racah or phase transforms. However, in what follows, in order to make exemplification simpler 
we shall limit our attention to rotation graphs only, i.e. 0„-graphs in which adjacent vertices differing only for a twist 
are identified, since they actually capture all the essential mathematical properties of our model. Fig.^ exhibits for 
the case n = 3 the local reduction of 03 when such identification is implemented (cfr. [0). 

Fig. 4 

Accordingly, the bonds of the rotation graph (that we still denote correspond to Racah transforms, possibly 
apart from weight/phase factors. Figj| shows an example of such reduced computational space for n = 3. 

Fig. 5 

The combinatorial structure of ©„ is fully determined by the identities connecting 6j symbols ||lo|] : 
i) the Bicdenharn-EUiot identity generates pentagon plaquettes in ©„: 




ii) Racah's identities generate triangles: 



here the spin variables {a, 6, c, . . . , x} ranges over all possible values in {O, |, 1, |, . . . } which obey the required 
triangular conditions, and R = a + b + c + d-\-e-\-f+p + q + r. We argue that the greater computational power of a 
quantum computer can be ascribed to the feature that its state space 'volume' grows very rapidly. Indeed the order 
(# of vertices) of ©„ is |©„| = {2n — 1)!! (~ n" for large n), whereas the diameter c)„ of C5„ grows approximately 
as nlnn ^ ln|C5„| JT^ . In the present scheme c)„ is an upper bound for the time-length (number of steps) of the 
computations machine 9Jt can perform (notice that for the full twist-rotation graph the cardinality is a factor 2" 
larger) . 
Universality 

Universality of dJl, being its computational space 25„, is guaranteed by Biedenharn-Louck's theorem (which plays 
the role of a sort of generalized Cayley's theorem): any unitary transformation corresponding to an operation of DJl 
can be reconducted to a finite sequence of operations in 0„. 

This gives an answer to the question raised by Feynman about universality j^] , explicitly defining the class of 'exact 
imitators' of any finite, discrete quantum system, with no need of resorting to the notions coming from the (inherently 
classical) Boolean circuit theory. 
Identification with Feynman's Q-Simulator 

VJl has all the requisites of the 'quantum simulator' as defined by Feynman Q: locality, refiected in the bracketing 
structure, which bears on the existence of local interactions; discreteness, both of the computational space and of 
'time', and universality. The time lapse from \in) to \out), as required in Feynman's scheme and as we shall see in the 
next section, is simulated both through the ordering induced by the graph combinatorial structure and by the number 
of computational steps; in other words, it is not only discrete but intrinsically inherent to the simulator structure [ p^ . 
It is the interplay between such space and time discreteness which gives rise in a natural way to entanglement, due 
to the clustering proper to the (non associative) Hilbert space tensor product structure generated by the recoupling. 



4 



DYNAMICS AND PROGRAMMING 

The above ingredients completely define the kinematical structure of 9Jl. Further notions are necessary to equip it 
with the dynamical structure necessary to make it operate. 

As in a classical Turing Machine computation is a map from the input data to the output state: U-p : \in) 
\out), where now the machine states are coded in vectors of the Hilbert spaces Ti^^ corresponding to the vertices of 
©„ and U-p is the class of unitary tranformations induced by the program V and defined by the corresponding 3nj 
symbols. 

The structure of computation in 9Jl is a generalization of the conventional Boolean scheme. To begin with, the 
coding language is based on an 'alphabet' consisting in all the (possibly different) values je , ^ = 1, . . . , rt + 1, of the 
coupled angular momenta, the intermediate variables fci, . . . , fc„_i, as well as the bracketing structure, and is therefore 
much more powerful and flexible. 

The program 7^(21) to perform algorithm 21 is the specification of a suitably designed {i.e. depending on 2t) ordered 
sequence of 'local alterations' of the alphabet elements in the running state, which play the role of gates. Such 
alterations are transforms of type (^), possibly accompanied by local permutations of labels and/or subtrees (moves 
of type (^, phase swaps) in 0„. We shall denote any such sequence by \ba), where index a, which keeps track 
of the given ordering, is such that Iba+i) is connected to \ba) by the elementary move corresponding to the local 
operation required by 21, while |bo} = \in) and Ifa^) = \out). L = L(7'(2l)) is the number of elementary steps required 
by program 7^(21) to complete algorithm 21. The associated lexicographically ordered sequence {\ba) \ a — 0, . . . , L} 
defines a directed path in 0„ of length L in one-to-one correspondence with the duration of V in units of its intrinsic 
discrete time step t. 

The associated unitary U-p(^^-j 

L-l 

{out\U-p^<ii^\in) = : '[[{ba+i\l^v{9i)\i'a) ■ , (4) 

Q=0 

where : • : denotes 'ordered product', is a sort of superselection rule which induces destructive interference of the 
forbidden {i.e. not leading to the correct result) paths in &„. Moreover, each elementary transfer matrix in (^ can 
be associated with a local hamiltonian operator 

{ba+i\l^v{<^)\ba) ^ exp{iH{ba,ba+i)T} . (5) 

It is worth noticing that 'local' is here intended with respect to the computational space ©„ of 971. When (^) is inserted 
in (^, hamiltonians H{ba, ba+i) generally do not commute with each other and are 'virtual', in the sense that they 
are generated by the machine dynamics in the course of the computation process. In the physical interpretation, 
however, they correspond to complex polylocal, many angular momentum binary interactions and simulate even more 
complex quantum physical systems (e.g., sets of interacting entangled fermions and/or bosons). 
Optimal computation and complexity 

Given two generic states |m), \out) in ©„, one can consider the Inverse Problem, namely read from the minimum 
length path the optimal algorithm that computes \out) as result of the application of some lA-p^t^^ to \in). The problem 
of finding the minimum-length path between two given vertices in 0„ is a 'hard' combinatorial problem, conjectured 
to be an NP-c problem (the question however is still open, iQ and references therein): we argue that the spin 
network simulator 971 may support algorithms to solve in polynomial time such kind of problems, because it is known 
[^5| that, at least in the case of unlabelled terminal nodes, the maximum distance between any pair of binary trees 
with N internal nodes is at most linear in N . 

Path-sum Interpretation and Topological Quantum Computation 

The dynamical behaviour of the spin network simulator is closely related to topological quantum field theories [ p^ . 
In particular, the sum over all ordered paths in between \in) and \out) of the functionals introduced in (^) has the 
form of a path- integral for a discrete topological quantum theory in (0-1-1) space-time dimensions. The topological 
(combinatorial) invariance is ensured by the equivalence of paths in 0„ under the set of identities for the 6j symbols 
discussed in Sect. II. Discretized models based on the recoupling theory of angular momenta have been extensively 
studied also in 3- and 4-dimensional quantum gravity on the grounds of the seminal paper [ pTf (see, e.g. and 
references therein). 

On the other hand, continuous gravity in (2+1) dimensions is well described by a Chcrn-Simons-Witten topological 
quantum field theory, whose basic objects are closed surfaces S of genus g. Freedman et al. in |^ resort just to 
the latter theory, by considering 'unitary topological modular functors' h, i.e. operations which - assigned a finite 
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dimensional Hilbert space to any such surface - connect diffeomorphic pairs S, S'. To each h there corresponds 

a transformation 7i(S) — > 7i(S'), reahzed in a quantum algebra. For special choices of this algebra, any such 
transformation is shown to be in one-to-one correspondence with a product of u elementary gate-operations, with Ti. 
interpreted as the computational space of a Quantum Circuit Model. The interger v, which measures the complexity of 
the corresponding 'computation', is linear in the length A of /i as a word in the standard generators (Dehn's twists) of 
the 'mapping class group' of S (whose compositions are cobordisms). Since A is in turn linear in the genus g, Freedman 
concludes that there exists a quantum circuit model that can efficiently simulate any topological modular functor in 
the given class. On the other hand, as the presentation of a 3-manifold by surgery and triangulations are equivalent 
p9[ , the approach described above can be in principle reconducted to a subclass of spin network simulators. This, 
roughly speaking, bears on the property that the cobordisms on a continuous manifold characterizing the modular 
functor can be translated into combinatorial moves over triangulations. 



CONCLUSIONS AND FURTHER DEVELOPMENTS 



We have exploited the kinematics and dynamics of a universal quantum simulator which encodes information in the 
full structure of the binary coupling schemes of angular momenta and manipulates such information by recoupling. 
The proposed model naturally exhibits all the characteristic features of the conventional model for quantum compu- 
tation, such as entanglement, intrinsic parallelism, tensor product structure of the state space, as well as the set of 
requirements proposed by Feynman as essential for the correct description of a quantum simulator. 

The model raises a number of intriguing questions, which of course demand extensive consideration; we mention 
just a couple of examples: 

Semiclassical Limit of the Simulator 

Given n, if the spin variables ji, . . . , jn+i, j, together with the intermediate fci, . . . , fc„_i, are all ^ 1 (in H units), 
each 6j symbol representing a Racah transformation can be approximated according to the asymptotics established 
in fill 




ji j2 ^1 \ ^ / ii h J5 1 / 1 

where V{T) is the euclidean volume of the tetrahedron T spanned by the six edges of 'length' £r — jr + \ i and 6r is 
the angle between the outward normals to the faces which share ^r (in the classical context canonically conjugate to 

This opens the intriguing possibility of bridging the quantum Turing Machine model to a classical counterpart and 
hence of interpreting in terms of classical actions the algorithms that it can solve. 
Mixed States Computation 

In principle, the whole conceptual scheme above can be reformulated in terms of density matrix formalism (i.e. 
resorting no longer to sharp eigenstates of the j^s but rather to generalized multipole moments |po| ). For example, 
in the simple case of two systems characterized by quantum numbers {(Ji), ji,mi , i = 1,2, where (ci) denotes all 
quantum numbers that are distinct from angular momenta eigenvalues, consider the Wigner coupling \{ai)jimi) 
|((72)j2W2) — > \io'iO'2jij2)j'm)- For given expansion of each single density matrix p^, i = 1,2, 

{{<^'t)jlm- \ p^\{ai)j^mi) = ^{(j[j[\\ Pt\c'iii)^^^C'^^l^^'^.^ , 

in terms of the reduced matrix elements and of the usual Clebsch-Gordan coefficients, one gets the expansion of the 
tensor product density matrix for the coupled system 



{{<^W2jlf2)fm'\pi P2\i<Jl<J2jlj2)jm) = ^{crW2jlf2f\\Pl ® P2\\<Jl<J2jlj2j) C] 



■jkj' 
mK,m' 

kK 



By standard methods of tensor operator theory pQ l , the reduced matrix elements above read 

f ^ ] k k 

K'J2j'i3'23'\\Pi®P2hi<J23ij2j)= Yl ^oky.o'A ki k2 k Wa'iiiipiiaiji)^; (a^j^||p2lk2j2),'^C^;^^^,« 

fclfc2KlK2 I j[ 32 j' I 
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where Wjkj'j' — [(2j + l)(2fc + l){2j[ + l)(2j2 + 1)] ^ , and the recoupling coefficients entering are 9j symbols. This 
comphcated expression gives the most general formula needed to describe any desired coupling by a sequence of binary 
couplings of density matrices, as required e.g. in a realistic quantum circuit implementation. 

We expect that the framework provided by the above remarks may permit including the environment in the picture, 
for example describing the simulator by a systems of pure angular momenta and the environment coupled one to 
another, either in terms of a density matrix or in the semiclassical approximation. 
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FIG. 1: Different rooted binary trees on (n + 1) =4 terminal nodes are depicted. Their (2n + 1) = 7 nodes are la- 
belled by angular momentum eigenvalues. The tree on the left corresponds to the binary bracketing \ji,j2,j3,j4 ; fci,fc2 ; j,rn) 
— * {{■'h{.'i2:i:i) f.^j4) J- The tree on the right corresponds to the binary bracketing \ji,j2,j3,j4',k'i,k'2;j,m) — > 

((jlJ2)^ (J3J4)^,)^.. 
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FIG. 2: The rotation operation on a portion of a labelled binary tree. The explicit expression of the Racah transform (g) 
relating the states associated with the trees depicted here reads \{a(bc)e)f ;m >— f^_i'^a+b+c+f ^^2fi ^ l)(2e + l)]^'^'^ 
{ c / e } I ( c)f]m> where the unitary matrix {6j} is the Racah- Wigner &j symbol of SU (2). 
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FIG. 3: The twist operation on a portion of a labelled binary tree. According to ^ the quantum state changes only by a phase 
factor. 
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FIG. 4: On the left there appears a local configuration representing six binary bracketings on four angular momentum variables 
{a, b, c, d}. Dotted lines represent twist operations (phase transforms on the corresponding state vectors), while the other edges 
are associated with rotations (Racah transforms between state vectors) . On the right the reduced configuration is shown, where 
now the bonds stands for one Racah transform plus some suitable phase/weight factors. The graph in Figj^ is built up taking 
into account this reduction procedure, and in particular its vertices {1,4,5} correspond to the vertices {Vi, V2,14} displayed 
here. 
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FIG. 5: The rotation graph (S3. Each vertex represents a binary coupling scheme of (n + 1) —4 angular momenta, two 
examples of which were given in Fig|l| The picture does not exhibit crossings but the vertices on the perimeter have to 
be identified through the action of the antipodal map (showing that (S3 is not planar). If we omit the intermediate labels 
of any binary bracketing of arguments {a,b,c,d}, the correspondences with the vertices are: 1 ^ (d(6(ac))); 2 ^ (6(d(ca))); 
3 ^ ((ac)(M)); 4 ^ (d(a(6c))); 5 ^ {d(c{ab))); 6 « (c(d(a6))); 7 ^ {{ab){cd)); 8 ^ (a(d(fec))); 9 ^ ((ad)(6c)); 10 « (a(6(cd))); 
11 ^ (a(c(M))); 12 ^ (c(a(M))); 13 ^ (c(b(ad))); 14 ^ (fe(c(da))); 15 ^ (6(a(cd))). 



